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The relativistic six-quark equations are found in the framework of the dispersion relation tech- 
nique. The approximate solutions of these equations using the method based on the extraction of 
leading singularities of the heavy hexaquark amplitudes are obtained. The relativistic six-quark 
amplitudes of dibaryons including the light quarks u, d and heavy quarks c, b are calculated. The 
poles of these amplitudes determine the masses of charmed and bottom dibaryons with the isospins 
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I. INTRODUCTION. 

Hadron spectroscopy has always played an important role in the revealing mechanisms underlying the dynamic of 
strong interactions. 

The heavy hadron containing a single heavy quark is particularly interesting. The light degrees of freedom (quarks 
and gluons) circle around the nearby static heavy quark. Such a system behaves as the QCD analog of familar 
_ hydrogen bound by the electromagnetic interaction. 
II ] The heavy quark expansion provides a systematic tool for heavy hadrons. When the heavy quark mass toq — >■ oo, 
f*) i- the angular momentum of the light degree of freedom is a good quantum number. Therefore, heavy hadrons form 
^ ■ doublets. For example, fib and fij will be degenerate in the heavy quark limit. Their mass splitting is caused by the 
chromomagnetic interaction at the order 0{l/mn), which can be taken into account systematically in the framework 

■ of heavy quark effective field theory (HQET) 

CsJ \ In 1977, JafFe Q studied the color-magnetic interaction of the one-gluon-exchange potential in the multiquark 
^ ' system and found that the most attractive channel is the flavor singlet with quark content v?(Ps^ . The same symmetry 
t ^ analysis of the chiral boson exchange potential leads to the similar result [SJ. 

The iJ-particle, ATJ-state and d\-fl may be strong interaction stable. Up to now, these three interesting candidates 
' of dibaryons are still not found or confirmed by experiments. It seems that one should go beyond these candidates 
• . and should search the possible candidates in a wider region, expecially the systems with heavy quarks, in terms of a 
1:0 ' more reliable model. 

There were a number of theoretical predictions by using various models: the quark cluster model 0,0] j the quark- 
delocation model [8, 9], the chiral SU(i) quark model [10], the flavor SU{i) skyrmion model [ll|. Lomon predicted a 
deuteronlike dibaryon resonance u sing R-matrix theory |12| . By employing the chiral SU{?>) quark model Zhang and 
Yu studied nn and SO states [H, [l4[ . 

In a series of papers [l5l - [T9j a method has been developed which is convenient for analyzing relativistic three- hadron 

■ systems. The physics of the three-hadron system can be described by means of a pair interaction between the particles. 
There are three isobar channels, each of which consists of a two-particle isobar and the third particle. The presence 
of the isobar representation together with the condition of unitarity in the pair energies and of analyticity leads to a 
system of integral equations in a single variable. Their solution makes it possible to describe the interaction of the 
produced par ticles in three-hadron systems. 

In Ref. [20| a representation of the Faddeev equation in the form of a dispersion relation in the pair energy in the 
two interacting particles was used. This was found to be convenient in order to obtain an approximate solution of the 
Faddeev equation by a method based on extraction of the leading singularities of the amplitude. With a rather crude 
approximation of the low-energy NN interaction a relatively good description of the form factor of tritium (helium-3) 
at low was obtained. 

In our papers [2l| - |23j relativistic generalization of the three-body Faddeev equations was obtained in the form of 
dispersion relations in the pair energy of two interacting quarks. The mass spectrum of S-wave baryons including u, d, 
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s quarks was calculated by a method based on isolating the leading singularities in the amplitude. We searched for the 
approximate solution of integral three-quark equations by taking into account two-particle and triangle singularities, 
all the weaker ones being neglected. If we considered such approximation, which corresponds to taking into account 
two-body and triangle singularities, and defined all the smooth functions of the subenergy variables (as compared 
with the singular part of the amplitude) in the middle point of the physical region of Dalitz-plot, then the problem 
was reduced to the one of solving a system of simple algebraic equations. 

In the present paper the relativistic six-quark equations are found in the framework of coupled-channel formalism. 
We use only planar diagrams; the other diagrams due to the rules of 1/Nc expansion [24i-,26j are neglected. 

The six-quark amplitudes of dibaryons are calculated. The poles of these amplitudes determine the masses of 
dibaryons. We calculated the contribution of six-quark subamplitudes to the hexaquark amplitudes. In Sec. II, 
we briefly discuss the relativistic Faddeev approach. The relativistic three-quark equations are constructed in the 
form of the dispersion relation over the two-body subenergy. The approximate solution of these equations using 
the method based on the extraction of leading singularities of the amplitude are obtained. We calculated the mass 
spectrum of 5'- wave bottom baryons with = i^, |^ (Tables HI HI]). In Sec. Ill, the six-quark amplitudes of 
hexaquarks are constructed. The dynamical mixing between the subamplitudes of dibaryons are considered. The 
relativistic six-quark equations are constructed in the form of the dispersion relation over the two-body subenergy. 
The approximate solutions of these equations using the method based on the extraction of leading singularities of 
the amplitude are obtained. Sec. IV is devoted to the calculation results for the dibaryon mass spectra and the 
contributions of subamplitudes to the hexaquark amplitude (Tables ITV l IV l IV H IVIip . In conclusion, the status of the 
considered model is discussed. 



II. BRIEF INTRODUCTION OF RELATIVISTIC FADDEEV EQUATIONS. 



In our papers, |21| - |23l |271 l28l | relativistic generalization of the three-body Faddeev equations was obtained in 
the form of dispersion relations in the pair energy of two interacting particles. The mass spectra of iS-wave baryons 
including u, d, s, c quarks were calculated by a method based on isolating of the leading singularities in the amplitude. 

We searched for the approximate solution of integral three-quark equations by taking into account two-particle and 
triangle singularities, the weaker ones being neglected. If we considered such an approximation, which corresponds 
to taking into account two-body and triangle singularities, and defined all the smooth functions at the middle point 
of the physical region of Dalitz-plot, then the problem was reduced to the one of solving a system of simple algebraic 
equations. 



In the paper 27| the relativistic three-particle amplitudes in the coupled-channels formalism are considered. We 



take into account the m, d, s, c, b quarks and construct the flavor-spin functions for the 35 baryons with the spin-parity 
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In the paper [28|, the relativistic equations were obtained and the mass spectrum of S'-wave charmed baryons was 
calculated. 

In the present paper, we will be able to use the similar method. In this case, we consider 35 baryons with the 
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spin-parity = ^ and Jp — ^ , which include one, two and three bottom quarks. We have considered the 23 
baryons with different masses (Tables HI HI]) . 

We calculate the masses of the bottom baryons in a relativistic approach using the dispersion relation technique. 
The relativistic three-quark integral equations are constructed in the form of the dispersion relations over the two-body 
subenergy. 

We use the graphical equations for the functions Aj{s, Sik)- In order to represent the amplitude Aj{s,Sik) in the 
form of dispersion relations, it is necessary to define the amplitudes of quark-quark interaction aj{sik). The pair 
amplitudes qq — ?> qq are calculated in the framework of the dispersion N / D method with the input four-fermion 
interaction with quantum numbers of the gluon. We use results of our relativistic quark model [29| and write down 
the pair quark amplitudes in the following form: 

aj{sik) = - r, 1 

1 - Bj{Sikj 

Aj(i,k) 
(mi+mfc)2 
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Here Gj is the vertex function of a diquark, which can be expressed in terms of the A^-function of the bootstrap 
N/D method as Gj = \^Nj, Bj{sik) is the Chew-Mandelstam function [s^, and p,j{sik) is the phase space for a 
diquark. Sik is the pair energy squared of diquark, the index Jp determines the spin-parity of diquark. The coefficients 
of Chew-Mandelstam function aj, /3j and 6j in Table Hill are given. Aj{i,k) is the pair energy cutoff. In the case 
under discussion the interacting pairs of quarks do not form bound states. Therefore, the integration in the dispersion 
integral ([2]) is carried out from (rrii + to Aj(i, k) {i, k — 1,2, 3). Including all possible rescatterings of each pair 
of quarks and grouping the terms according to the final states of the particles, we obtained the coupled systems of 

integral equations. For instance, for the EjJ" with Jp = the wave function is ip^+ = \ / ^{u ^ d ^ b 1} — y ^{u ^ d I 



b^} — y I d t 6 t}. Then the coupled system of equations has the following form: 

' Ai{s, S12) = A6i(si2)ii(si2) + A"i(si2) [3^16(5, S13) + jAQb{s, 813)+ 

+ i^l^(s,S23) + |^i.(s,S23)] 
^l''(s,Si3) = A6ib(si3)Lib(si3) +Kib{si3) [i^l(s,Sl2) - \Aib{s,Si2) + 

+ |^Of('S, S12) + 5^l(s,S23) - zAib{s,S23) + jAQb{s,S23)] 
Ao''{s,S23) = A ^Qb (523)^0" (S23) +Kot,{s23) [^Ai{s,Si2) + jAib{s, 812) + 
+ jAQi{s,Si2) + ^Ai{s,Si3) + jA^b{s,Si3) + jAQb{s,Si3)] ■ 

Here the function Lj{sik) has the form 



(4) 



The integral operator Kj(sik) is 



(5) 
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(mi+mt.)2 -1 



The function bj{sik) is the truncated function of Chew-Mandelstam: 



CX) 

bj{sik) = J 

z is the cosine of the angle between the relative momentum of particles i and k in the intermediate state and the 
momentum of particle j in the final state, taken in the cm. of the particles i and k. Let some current produces three 
quarks with the vertex constant A. This constant do not affect to the spectra mass of bottom baryons. By analogy 
with the state, we obtain the rescattering amplitudes of the three various quarks for the all bottom states. 
Let us extract two-particle singularities in Aj{s, Sik)'- 
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(7) 



. , N Oij{s,Sik)hj{sik)Gj{sik) , , 

A,J[s,S^k) = ^-7 T , 8) 

ctj{s,Sik) is the reduced amplitude. Accordingly to this, all integral equations can be rewritten using the reduced 
amplitudes. The function aj{s, Sik) is the smooth function of Sik as compared with the singular part of the amplitude. 
We do not extract the three-body singularities, because they are weaker than the two-particle singularities. For 
instance, one considers the first equation of system for the Y.^ with JP — 
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The connection between s'12 and 5^3 is [2l| 



, 2,2 {^'12 +ml~s) (s'i2 +ml- mf) 
3i3 = mi-fm3 — 



± 



2S]^2 



'-\/is'l2 - (^1 + ^2)^) (S'l2 - (™l - "12)2) {s[2 - (V^ + ^3)2) (s'la -iV^S- 7713)2) . (10) 



The formula for S23 is similar to (|10p with replaced by z — >■ ~z. Thus Aib{s, s[^) + Aib{s, S23) must be replaced by 
2Aib{s, s[^). Aj(i, k) is the cutoff at the large value of Sik, which determines the contribution from small distances. 

The construction of the approximate solution of coupled system-equations is based on the extraction of the leading 
singularities which are close to the region Sik = {mi + rnk)"^ [3l| . 

We consider the approximation, which corresponds to the single interaction of the all three particles (two-particle 
and triangle singularities) and neglecting all the weaker ones. 

The functions aj(s, Sik) are the smooth functions of Sik as compared with the singular part of the amplitude, hence 
it can be expanded in a series at the singulary point and only the first term of this series should be employed further. 
As So it is convenient to take the middle point of physical region of the Dalitz-plot in which z = 0. In this case, we 

get from Sik = so = , where uiik = "''t™'' . We define functions aj{s,s^k) and bj{sik) at the point 

So- Such a choice of point so allows us to replace integral equations (|3]) by the algebraic couple equations for the state 
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(11) 



ai (s, So) = A + i a^b (s, sq) (s, sq) + | agb (s, so) /lo" (s, so) 

aii,(s,So) = A + q;i(s,So)/ih(s,so) 

-iaib(s,so)/i!>i!.(s,so) + |aob(s,so)/i!.o''(s,so) \"^X°o) 
aob(s,so) = A + Q!i(s,so)/obi(s,so) 

+i aib(s,so)/obib(s,so) '^-^j + ^aoi'(s,so) /obib(s, sq) . 

The function Ij^j^{s,sq) takes into account singularity which corresponds to the simultaneous vanishing of all 
propagators in the triangle diagram. 
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The Gj{sik) functions have the smooth dependence from energy Sik [25|, therefore we suggest them as constants. 
The parameters of model: gj vertex constant and Aj cutoff parameter are chosen dimensionless; 



rrii+mk ^ . 4Aj 

9j 7^ ^Ji — 7 
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(13) 



Here rrii and nik are quark masses in the intermediate state of the quark loop. We calculate the coupled system 
of equations and can determine the mass values of the T,'^ state. We calculate a pole in s which corresponds to the 
bound state of the three quarks. 

By analogy with EjJ^-hyperon we obtain the systems of equations for the reduced amplitudes of all bottom baryons. 

The solutions of the coupled system of equations are considered as: 



a,/ 



Fj{s,Xj) 
D{s) 



(14) 



where the zeros of the D{s) determinate the masses of bound states of baryons. Fj{s, Aj) are the functions of s and 
Xj. The functions Fj{s, Aj) determine the contributions of subamplitudes to the baryon amplitude. 

In quark models, which describe rather well the masses and static properties of hadrons, the masses of the quarks 
usually have the similar values for the spectra of light and heavy baryons. However, this is achieved at the expense 
of some difference in the characteristics of the confinement potential. It should be borne in mind that for a fixed 
hadron mass the masses of the dressed quarks which enter into the composition of the hadron will become smaller 
when the slope of the confinement potential increases or its radius decreases. Therefore, conversely, we can change 
the masses of the dressed quarks when going from the spectrum of light baryons to the heavy baryons, while keeping 
the characteristics of the confinement potential unchanged. We can effectively take into account the contribution of 
the confinement potential in obtaining the spectrum of S'-wave heavy baryons. 

In the case of b quark we have used two new parameters: the cutoff of the bb diquark Af, — 5.4 and the coupling 
constant Of, = 1.03. These values have been determined by the 6-baryon masses: 1+ = 5.808 GeT^ and 3+ = 
5.829 Gey. In order to fix mt = 4.840 GeT^ we use the 6-baryon mass 5.829 GeV. We represent the masses of all S'- 
wave bottom baryons in the Tables HI HIl The calculated mass value M^^ 1+ = 5.624 GeV is equal to the experimental 

data (13, the mass value M~a 1+ = 5.761 GeV is close to the experimental one. We neglect with the mass distinction 

^sb 2 

of u and d quarks. The estimation of the theoretical error on the bottom baryon masses is 2 — 5 MeV. This result 
was obtained by the choice of model parameters. 

In our model the spin-averaged mass of the states and is predicted to lie around to 250 MeV above M^^ . 
The relativistic corrections are particularly important for the splitting between and f2{, baryons. 

In the context of SJ, and masses, it is worth mentioning two relations among bottom baryons which incorporate 
the effects of SU (3) / breaking: 
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[M (El) - M (Sfc)) + {M (Ql) - M (nt)) - 2 (M (S^) - M (4)) = , (15) 
M(5]b) + M(f7fc)-2Af(4) = 0. (16) 



The sign in our prediction is 



M (Efc*) - M (Efc) < M (ni) - M {ilb) . (17) 



This inequaUty is not predicted by other recent approaches [331 l34l | , which predict a fib sphtting smaller than a 
Eb sphtting. This suggests that the sign of the SU{3) symmetry breaking gives information about the form of the 
confinement potential. 

We have used the mb/mc = 2.95 in the Tables HI and lU (similar to the Ref. [ssjV 

III. SIX-QUARK AMPLITUDES OF THE HEXAQUARKS. 

We derive the relativistic six-quark equations in the framework of the dispersion relation technique. We use 
only planar diagrams; the other diagrams due to the rules of 1/Nc expansion [23 - l26| are neglected. The current 
generates a six-quark system. The correct equations for the amplitude are obtained by taking into account all 
possible subamplitudes. Then one should represent a six-particle amplitude as a sum of 15 subamplitudes: 



i<3 
1,3 = 1 

This defines the division of the diagrams into groups according to the certain pair interaction of particles. The 
total amplitude can be represented graphically as a sum of diagrams. We need to consider only one group of diagrams 
and the amplitude corresponding to them, for example Ai2. We shall consider the derivation of the relativistic 
generalization of the Faddeev-Yakubovsky approach. In our case, the low-lying dibaryons are considered. We take 
into account the pairwise interaction of all six quarks in the hexaquark. 

For instance, we consider the state AA;, with ^ = |, = 2+ and quark content uuuudb. The set of diagrams 
associated with the amplitude A12 can further be broken down into eight groups corresponding to subamplitudes: 

41"" /lO"'' /lO"'' AO'"' .ii'^o*"' 4l""0"'' 4r"'0'"' .i^^i^^Qdb 

The amplitude A\ (s, S12345, S1234, S123, S12) consists of the three color sub-structures: the diquark 1"" in the color 
state 3c, the quarks 3, 4 in the color state 3c x 3c = 3c -I- 6c, and the quarks 5, 6 in the color state 3c x 3c = 
3c -I- 6c. Then we consider the total color singlet: 3c x 3c x 3c = Ic + 8c + 8c + 10*. The dibaryon amplitude 
Al ° (s, S12345, S1234, S12, S34) contains the following sub-structures: the two diquark 1"" and 0"'' in the color state 
3c and the two quarks in the color state 3c. Then the dibaryon amplitude is the total color singlet. The amplitude 
2345, S12, S34, S56) consists of the three diquark structures in the color state 3c. Therefore the total 
color singlet can be constructed. For the others amplitudes color singlet also can be found. 

The system of graphical equations (see for example equation for the amplitude A2 ° for the state AAf, with I = ^ 
and = 2+ at the Fig. 1) is determined by the subamplitudes using the self-consistent method. The coefficients 
are determined by the permutation of quarks. 

In order to represent the subamplitudes Af"', Af\ Al''\ ^r'^o"-^^ ^^u^^ut ^ ^^„upd6^ ^^^u^^^^dt ^^^^ 

of a dispersion relation, it is necessary to define the amplitude of qq and qQ interactions. We use the results of our 
relativistic quark model (29j and write down the pair quark amplitudes in the form: 

an(Stk) = -. 5-7 — 7, (19) 



B„(s.fe)= / ^^liifkMlfkl, (20) 

J 1" S^fc - S,fc 

(mi+mfc)2 
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[nil + mfc)2 
\/(sifc - {rrii +mkY){sik - {rrii - rukY) 



(21) 



The coefficients a{n, J^'^), (3{n,J^'-^) and S{n,J^'~') are given in Table IVllIl Here n = 1 coresponds to qq and 
qQ-pairs with =0+, n — 2 corresponds to qq and gQ-pairs with = 1+. 

The coupled integral equations correspond to Fig. 1 can be described similar to 36]. Then we can go from the 
integration of the cosine of the angles dzi to the integration over the subenergies. 

Let us extract two- and three-particle singularities in the amplitudes A}"", A^" , A^" , , yli""o" ^ Aj""^" , 
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[1 - Biu.(si2)][l - B1...(S34)][1 - Bodb{s56)] 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 

(29) 



We used the classification of singularities, which was proposed in paper [31|. Using this classification, one defines 
the reduced amplitudes ai, 02, a^, as well as the _B-functions in the middle point of physical region of Dalitz-plot at 
the point so- 

Such choice of point sq allows us to replace integral equations (AA;,, I — ^, — 2+) by the algebraic equations 



af " = A + 4 af "/i (1™1™) + 2 af /i (1""0"'*) + 2 af h (1™'0"'') + 4 "°"'/2(l™l™0"'') 
+ 4a^""°"''/2(r"r"0"''), 

= A + 3af "/i(0"'^l"") + i h{Q'"^Q'"^) -f af"'/i(0"'^0"'') -f af''/i(0"''0'^^) 
6ar"°"'/2(0"''l™0"'') + 3ar"°'"'^2(0"'^r"0''^) , 



(30) 



(31) 



af = A + 3af "7i(0"''l"") + /i(0"^0"'*) + Saf /i(0"''0'''') + /i(0"^0'''') 

+ 6ar"°"'"/2(0"''r"0"'') + 3af"°"'/2(0"''l""0'"'), 

= A + 4ar7i(0''^0"<^) + 4ar7i(0*0"^), 

'"0"' = A + af" (2 /3(i««o"'^l"") + 2 /4(l""0'"'l"")) + 2 af" J3(l""0"''0"'') + af 74(0"''l""0"'') 

+ af 74(0"'*1""0'^^) + 2 a^""°"'77(0'"^l""0'"^l"") + ai"""'" (2 /5(i™o"<^l""0"'') 

+ 2 76(l""0"''l""0"*)) + a^""°'"'(75(0"''l""l""0''^) + 2 76(1""0"'^1""0''^) + 2 77(l""0"''l""0'^'')) 

+ 2 4""i""°'"'78(l""0"''l""l"V') , 

.UQUb ^ ^ _^ ^1- ^2 73(i«»o"*l"") + 2 74(l""0"''l"")) + 74(0"*1""0"'^) + 2 af 73(1""0"''0"*) 

+ 74(0"''1''"0'^'') + a^""°"''(2 75(l""0"''l""0"'^) + 2 76(1""0"''1""0"'')) 

+ 2ar"°"''77(0"''l""0"''l"'')) + ar"°'''(75(0"^l""l""0'^'') + 2 76(r"0"''l""0''^) 

+ 2 77(1""0"^1""0'*'')) + 2 aP^""°'"'78(l""0"^l""l""0'*'') , 

""o"' = A + 4af "74(l""0'^''l"") + (2 73(1"V0"'') + 2 74(0'^''1""0"'')) + af (2 73(r"0''''0"'') 

+ 2 74(0''^1""0"^)) + ar''°"''(4 76(l""0''^l""0"'') + 477(0''^l""0"''l"")) 

+ a^""°"''(476(l""0''''l«"0«'') +4 77(0'^''l"«0"''l"")) , 

""0" = A + 4ap79(l"''l"V''l"") +4af'' 79(1"V*1""0'"') +40?"' 79(1""0''^1™0"'') 

+ 8Q!r"°"''7io(l""l""0''^l""0"'') + 8Q!r"°"''7io(l""l""0''''l""0"'') , 



where A, are the current constants. We used the functions h, I2, I3, h, h, h, I7, h, h, ho' 
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-1 -1 



X 




hiijk) = 



(mi+m2)^ 



(m3+m4)^A- 



+rr " 

h{ijk) = h{ik) 



23 

(m3+m4)^ —1 —1 



h{ijkl) = hiikl), 
h{ijkl) = h{ik)-h{jl), 



f ds'^^G]{sf)pj{s'3^) 1 JdziiJ) J' dz2{7) JdzaiJ) 
J TT s's^-S3i 2ttJ 2 J 2 J 2 

(ma+mi)^ -1 -1 -1 

Z4(7) + 

X j dzi{7) 



VI - ^?(7) - ^1(7) - zl{7) + 2.1(7)^3(7)^4(7) 

24(7)- 



Is{ijklm) = 



1 - Bfe(sy 1 - ' 

(mi+m2)^A. 



i?,(sl2)S,(s34) J ;r S[,-Sl^ 

(mi+m2)^ 

(m3 + m4)^Aj- 

; rf5^4G,'(«o')ft(44) 



o' _ o34 

S34 



(m3+m4)^ 

+ 1 +1 +1 24(8)+ +1 ze(S) + 

-1 -1 -1 24(8)- -1 26(8)- 

1 

- zt{8) - zi{8) - z|(8) + 2zi (8)^3(8)24(8) 
1 

^1 - z|(8) - z|(8) - .2(8) + 2z2(8)z5(8)z6(8) 
111 

1 - Bfe(s'l5) 1 - BK43) 1 - 5m(si6) ' 

hiijkl) = h{ijl), 
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Iio(ijklm) 



(mi+n.2)^Aj 



'12 ^0 



X 

'34 ~ ^0 " ''56 ~ ''0 

+ 1 +1 +1 +1 Z5(10) + 

1 J dzi(lO) J dz2{lO) J dz3(10) J rfz4(10) J ^^^^^^^ 

-1 -1 -1 -1 25(1-)- 

1 

Vl - z?(10) - ^I(IO) - z|(10) + 2zi(10)z4(10)z5(10) 
1 1 



l-i?K.s'23)l-S™K5) 



(47) 



where i, j, k, I, m correspond to the diquarks with the spin-parity = 0+, 1^. 
The solutions of the system of equations are considered as: 

a,(s) =F,(s,A,)/i^(s), (48) 
where zeros of D{s) determinants define the masses of bound states of dibaryons. 

IV. CALCULATION RESULTS. 

The model in question has three parameters of previous model (36| : gluon coupling constants go — 0.653 (diquark 
0''") and gi = 0.292 (diquark 1+), cutoff parameter A — 11. We used to the cutoff A^h — 4.43 and the cutoff A^c = 5.18, 
which are determined by M = 7300 MeV^ (threshold 7315 Me^) and M = 4100 AfeF (threshold 4130 Afe^). The 
experimental data is absent, therefore we use the dimensionless parameters, which are similar to the Eq. It 
allows us to calculate the mass spectra of qqqqqQ states. 

The quark masses of the model are = 495 M eV, rric = 1655 MeV, mi, = 4840 MeV . The estimation of theoretical 
error on the S'-wave hexaquarks masses is 1 M eV. This results was obtained by the choice of model parameters. 

We have considered 38 dibaryons with content qqqqqQ, q = u,d, Q = c,b. The masses of dibaryons with I — ^, 
|, I and spin-parity = 0+, 1+, 2+ in the Tables JV] and |V] are given. The lowest mass for the qqqqqb states is 
M = 5700 MeV. The lowest mass for the qqqqqc states is M = 3475 MeV. 

The relativistic six-body approach possesses the dynamical mixing and allows us to calculate the contributions of 
the subamplitudes to the hexaquark amplitude (Tables IVll I VII|) . The calculated dibaryon subamplitudes A2 present 
the main contributions to the hexaquark amplitude (about 80 percents). 

In a strongly bound systems, which include the light quarks, where p/m ^ 1, the approximation by nonrelativistic 
kinematics and dynamics is not justified. 

In our paper, the relativistic description of three-particles amplitudes of bottom baryons are considered. We take 
into account the u, d, c, b quarks. The mass spectrum of S'-wave bottom baryons with one, two, and three b quarks 
is considered. We use only two new parameters for the calculation of 23 baryon masses. The other model parameters 
in the our papers [2l| - |23| are given. 

We have predicted the masses of baryons containing b quarks using the coupled-channel formalism. We believe that 
the prediction for the S'-wave bottom baryons based on the relativistic kinematics and dynamics allows as to take 
into account the relativistic corrections. In our consideration, the bottom baryon masses are heavier than the masses 
in the other quark models [13, HI] . 

Our model is confined to the quark-antiquark pair production on account of the phase space restriction. Here ruq is 
the "mass" of the constituent quark. Therefore the production of new quark-antiquark pair is absent for the low- lying 
hadrons. 

We manage with quark as with real particles. However, in the soft region, the quark diagrams should be treated 
as spectral integral over quark masses with the spectral density p{m?): the integration over quark masses in the 
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amplitudes puts away the quark singularities and introduces the hadron ones. One can believe that the approximation: 
p[m?) — ^ 5{m? — m^) could be possible for the low-lying hadrons. 

The authors of Ref. [s^ calculated the energies of the baryon-baryon threshold as a function of the flavor-symmetry 
breaking parameter 5 = \ — . The binding energy is obtained and the possible candidates for stability under strong 
interactions is commented on. In the recent papers (sol . |40| , developed for the baryon-baryon interactions in lattice 
QCD, the flavor-singlet H dibaryon is studied. The results of the lattice QCD calculations presented the first clear 
evidence for a bound state of two baryon directly from QCD. 
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TABLE I: Bottom baryon masses of multiplet Parameters of model: quark masses mu,d = 495 MeV, me = 770 MeV, 
rric = 1655 MeV, mt = 4840 A/eV; cutoff parameters: Xq = 10.7 (g = u,d,s), Ac = 6.5, At = 5.4; gluon coupling constants: 
50 = 0.70, gi = 0.55 with J" = 0+ and 1+, ffc = 0.857, gt = 1.03. 



Baryon 


Mass (GeV) 


Mass (GeV) (exp.) 


Eft 


5.808 


5.808 




5.624 


5.624 




5.761 


5.774, 5.793 


' — 'sb 


6.007 




ssb 


6.120 






6.789 




^cb 


6.818 






6.798 






6.836 




^ccb 


7.943 




S(,6 


10.045 






9.999 




iicbb 


11.089 





TABLE II: Bottom baryon masses of multiplet Parameters of model: quark masses mu,d = 495 MeV, rus = 770 MeV, 
rric = 1655 MeV, mb = 4840 Mel^; cutoff parameters: \q = 10.7 {q = u,d,s), Ac = 6.5, Xb = 5.4; gluon coupling constants: 
go = 0.70, gi = 0.55 with = 0+ and 1+, 5c = 0.857, gb = 1.03. 



Baryon 


Mass (GeV) 


Mass (GeV) (exp.) 




5.829 


5.829 




().()()() 




^ssb 


6.220 




Hcb 


6.863 




^scb 


6.914 




^ccb 


7.973 






10.104 




i^sbb 


10.126 




iicbb 


11.123 




0,bbb 


14.197 





TABLE III: Coefficients of Ghew-Mandelstam functions. 
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TABLE IV: S-wave charmed dibaryon masses. Parameters of model: cutoff A = 11.0 and Aqc = 5.18, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses rUq = 495 MeV and rUc = 1655 MeV. 



T 


7 

J 


Dibaryons (quark content) 


iViarSS i^iviev ) 


5 
2 


n 
U 


AS* (uuu uuc) 


41UU 




1 o 

i, 2 


AV AV* 

ZAZjc, ^-^c (uuu uuc) 


/1 1 nn 
4iUU 


■A 
2 


u 


AS* (^uuu udc -\- uud uuc) 


oo / u 






iVzjc (uwa uucj 


ooyy 




i 


Azjc, Azjc (unn uac-\-uua uuc) 


oo/U 






AAc {uuu udc) 


oyzu 






ly i^c-, iVzjc (iiM tfucj 


ooyy 




2 


Azjc, AZjc (mmm udc + uud uuc) 


3746 






AAc {uuu udc) 


oyzu 






iVZjo (ttttCl MMCj 




1 
2 







3475 






1\ I\c \UUCl UdC) 


oozy 






NYlc {uud udc + itcZcZ uuc) 


o4oU 




1 


AEc, AE* (uuu ddc + uud udc + udd uuc) 


3475 






TV" Ac, A Ac (uud udc) 


3629 






TVEc, NY^c {uud udc + udd uuc) 


3480 




2 


AEc, (uuu ddc + uud udc + udd uuc) 


3863 






AAc (uud udc) 


3937 






A/"E* (uud udc + udd uuc) 


3870 
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TABLE V: S-wave bottom dibaryon masses. Parameters of model: cutoff A = 11.0 and Aqb = 4.43, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses niq = 495 MeV and nib = 4840 MeV. 



T 
1 


7 

J 


Dibaryons (quark content) 


Mass (^iviev ) 


5 
2 


u 


AS^ ('U'U'U 


/ ouu 




i, z 


AZjb, (ttWW UUOj 


/oUU 


3 
2 


U 


AS^ {uuu udb + nnd uitfe) 


oyoo 






iVzjb (liua uuo) 


DZOO 




i 


Aijb, [UUU uao + tiwa nnoj 


oyoo 






AA^ {uuu udb) 


DyZD 






lyi^bt jVzj(, (^'U'ua titiOj 


OZOO 




1 


A2j(,, ZaZj;, (MMM Mao + WUa MMOj 


o4dU 






AA^ (uuu udb) 


oyzD 






/"^.^.J 
iVZji, (MMffl MUOj 


DoUU 


1 

2 


U 


AE^ (mmm dd6 + ititd ud6 + udd uub) 


KTAA 
(UU 






iVAfo yuuOb uao) 


D14Z 






NJ^b {uud udb + udd uub) 


/ zo 




1 


AEfe, AS^ (www ddb + uud udb + udd uub) 


5700 






NAb, AAb {uud udb) 


6142 






NTfb, NT.I {uud udb + udd uub) 


5723 




2 


AEj,, AEj (www ddb + uud udb + udd uub) 


6744 






AAb {uud udb) 


6928 






N'El {uud udb + udd uub) 


6755 



TABLE VI: AEc, AE* {3475 MeV) {IJ = §1). Parameters of model: cutoff A = 11.0 and Age = 5.18, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses iriq = 495 MeV and rric = 1655 MeV . 



Subamplitudes 


Contributions, percent 


Ar 

Al"" 


4.98 
4.20 


Af' 


3.81 


Ar 


0.82 


Af 


0.23 


AC''' 


9.58 




17.73 




1.17 




1.69 


^2 


1.43 




46.54 




3.23 


Ql.<iQ<iC 

^2 


3.11 


^3 


0.43 




1.06 




14.05 




84.46 




1.49 
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TABLE VII: AE6, AS^ (5700 MeF) (7J = ^l). Parameters of model: cutoff A = 11.0 and A,;, = 4.43, gluon coupling 
constants go = 0.653 and gi = 0.292. Quark masses ruq = 495 MeV and mt — 4840 MeV. 



Subamplitudes 


Contributions, percent 




5.79 
4.89 


Af" 


2.67 


Al 
At 

Ar'"" 


0.11 
0.02 
11.04 


Ar'"' 


19.82 


^2 


0.13 
0.21 
U.iD 


A2 


04. it) 


A2 

^2 
^3 
^3 


0.39 

0.06 
0.14 




13.47 




86.33 


E^3 


0.20 



TABLE VIII: Vertex functions and Chew-Mandelstam coefficients. 



i 






Pi 


5i 


0+ 


3 (33fci) 


1 
2 


1 ("ifc-m,)^ 





1+ 


2g 


1 


4mfcmi 1 


1 (wifc-mi)^ 


3 


3 


3(mfc+m,()^ 6 


6 (mfc+m;)^ 



